Role of the phase-matching condition in non-degenerate four-wave mixing in hot 
vapors for the generation of squeezed states of light 
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Non-degenerate forward four-wave mixing in hot atomic vapors has been shown to produce strong 
quantum correlations between twin beams of light [McCormick et al, Opt. Lett. 32, 178 (2007)], 
in a configuration which minimizes losses by absorption. In this paper, we look at the role of the 
phase-matching condition in the trade-off that occurs between the efficiency of the nonlinear process 
and the absorption of the twin beams. To this effect, we develop a semi-classical model by deriving 
the atomic susceptibilities in the relevant double-lambda configuration and by solving the classical 
propagation of the twin-beam fields for parameters close to those found in typical experiments. These 
theoretical results are confirmed by a simple experimental study of the nonlinear gain experienced 
by the twin beams as a function of the phase mismatch. The model shows that the amount of 
phase mismatch is key to the realization of the physical conditions in which the absorption of the 
twin beams is minimized while the cross-coupling between the twin beams is maintained at the 
level required for the generation of strong quantum correlations. The optimum is reached when the 
four-wave mixing process is not fully phase matched. 

PACS numbers: 42.50.Gy,42.65.Yj,42.65.Hw 



I. INTRODUCTION 

Continuous-variable entanglement can be generated 
deterministically with a phase-insensitive optical ampli- 
fier. For a gain larger than 1, the system produces a 
two-mode squeezed state where the signal and the idler 
(here referred to as probe and conjugate respectively) 
display EPR-type entanglement pQ. Such an amplifier 
can be realized using a nonlinear optical process such 
as parametric down-conversion 2. 3 or four- wave mix- 
ing (4WM) In real physical systems, the presence of 
absorption reduces the amount of quantum correlations 
which can be generated, and although 4WM in atomic 
vapors can lead to large gains, resonant atomic processes 
are also responsible for losses, which limit the amount of 
observable squeezing. Recently, a configuration in 4WM 
was found which reduces absorption 5, 6 , and generates 
large degrees of squeezing [7HTU], 

There have been a number of reasons put forward to 
explain this success. The main one is the nature of the 
nonlinearly, which is based on coherence effects between 
the hyperfine electronic ground states rather than on the 
saturation of a transition of a two- level atom ITT] . 
Indeed, avoiding a large atomic population in the ex- 
cited state is key to the reduction of the noise associ- 
ated with spontaneous emission. More specifically, it was 
pointed out that the Dl line of alkali atoms is particu- 
larly amenable to the establishment of a ground state 
coherence [12] , We show here that the production of 
squeezing is also due in great part to a judicious choice 
of the parameters that most greatly influence the phase- 
matching condition of the nonlinear process, specifically 
the relative frequencies of the beams and the angle be- 



tween the beams. Maybe surprisingly, the highest levels 
of squeezing are achieved when the system is not fully 
phase matched. 

The paper is divided as follows. In section [TTJ we the- 
oretically study 4WM in an atomic vapor in a double- 
lambda configuration where both pumps are detuned 
from the atomic resonance. From the atomic suscep- 
tibilities we evaluate the impact of the phase-matching 
condition on the gain and the absorption in the forward- 
4WM configuration. In section |III| we report on a sys- 
tematic experimental study of the phase-matching con- 
dition which confirms the findings of section [n] In sec- 
tion |IV| we extend our model to take into account the 
Doppler effect due to the thermal motion of the atoms. 
Finally in section [V] we discuss the impact of our the- 
oretical and experimental findings on the possibility of 
generating strong two-mode squeezing with 4WM in a 
hot vapor. From the model we deduce the best param- 
eters in terms of beam geometry and beam detunings, 
and compare them to the recent squeezing experiments 
of Refs. EHini 



II. THEORETICAL MODEL OF 
NON-DEGENERATE FOUR- WAVE MIXING 

When describing nonlinear media in the presence of 
off-resonant fields, it is common to separate the response 
of the system into a linear part and a nonlinear part [13] . 
The linear contribution leads to an index of refraction 
which modifies the linear dispersion relation for each of 
the individual light fields in the medium. The nonlinear 
part acts as a perturbative source term in the propaga- 
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tion equations. It enables energy transfer between light 
fields for those configurations where the phase-matching 
condition is fulfilled, that is to say when the total wave- 
vector of the waves giving up energy equals the total 
wave- vector of the waves receiving the energy. The rel- 
evant wave-vectors are those in the medium. They are 
equal to the wave-vectors in vacuum times their corre- 
sponding indices of refraction. 

In resonant media, such as atomic vapors excited close 
to atomic transitions, the medium is strongly perturbed 
by the presence of the light. For instance, optical pump- 
ing by a strong pump beam can affect the atomic popu- 
lations in the hyperfine levels, leading to a strong change 
in the index of refraction seen by a weaker beam. In 
this case, the usual expansion separating the linear and 
the nonlinear responses may not be appropriate. Instead 
we consider an expansion of the nonlinear polarization to 
first order in the electric fields of the weak beams and to 
all orders in the electric fields of the pump beams [6]. 
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Figure 1. Double-lambda scheme on the Dl line of 85 Rb. 
The hyperfine splitting of the excited state is not resolved 
due to Doppler broadening. The lambda transitions, detuned 
by Ai and A2, are driven by pump fields with resonant Rabi 
frequencies fii and Q2, respectively. Note that the pumps can 
actually be a single laser beam. 

We consider the double-lambda configuration, shown 
in Fig. [l] which was used to demonstrate intensity- 
difference squeezing and quadrature entanglement in a 
vapor of 85 Rb [T2J [15]. A non- degenerate 4WM para- 
metric process drives an atom from one of the hyperfine 
ground states to the other hyperfine ground state and 
back to the initial state. In the process, two pump pho- 
tons are converted into two twin photons, called probe 
and conjugate, with wave-vectors k p and k c in vacuum, 
and frequencies uj p and uj c . The non-linearity originates 
in a strong coupling between the probe and the conjugate 
fields mediated by the coherence of the electronic ground 
states Following the usual experimental configura- 
tion, we further assume that the two pump photons come 
from a single field, of wave-vector k in vacuum and fre- 
quency (jJq- We denote 5 the two-photon detuning of the 
pump and the probe: S = ujq — uj p — whf where whf is 
the hyperfine splitting of the ground state. The natural 
linewidth of the excited state is 7 — 2ir x 6 MHz. In the 
rest of the paper, we vary 5 by changing the frequency of 



the probe. The 4WM resonance occurs roughly at <5 = 
(more on this below). 

For simplicity, we assume that the coupling strengths 
of the pump to both transitions are equal, correspond- 
ing to a resonant Rabi frequency O = £7i = £^2- Fol- 
lowing Ref. we calculate the atomic susceptibilities at 
the probe and conjugate frequencies in the limit of weak 
probe and conjugate fields. The detailed calculation is 
developed in the appendix. The susceptibilities are de- 
rived by calculating the steady-state value of the density 
matrix of a 4-level system interacting with the four fields 
of the double lambda. The atomic polarization at the fre- 
quencies of the probe and the conjugate is proportional 
to the average oscillating atomic electric dipole at those 
frequencies and therefore to the off-diagonal components 
of the density matrix corresponding to these transitions. 
To first order in the probe and conjugate fields, the 
atomic polarization is described by two direct suscep- 
tibilities, Xpp an d Xcc, and two cross-suscep tibilities, Xpc 
and Xcp = Xpci given by Eqs. (A.13-A.16). The cross- 
susceptibilities are responsible for the 4WM. 

The propagation equations for the slowly varying en- 
velopes of the probe and conjugate fields £ v an d £ c , using 
the polarization expressions (A. 11 ) and (A. 12), are given 
in steady state by: 



dz tp ~ 2e 1 p) 

& c ike r>l \ -ik c 

dz £c= 2V P{uJc)e 



(1) 
(2) 



Furthermore, if we consider the case of co-propagating, 
or nearly co-propagating, beams along the z axis, these 
equations read: 

§~/p = ^Xpp^pK + ^Xpc^pY^E: (3) 

§-/c = ^xceMfc + ^ Xcp KK Afc -£;, (4) 

where Ak z is the projection of the geometric phase- 
mismatch Ak = 2ko — k p — k c on the z axis, and the con- 
servation of energy imposes the condition uj p + cj c = 2w - 
In the low pump depletion limit, which is usually ex- 
perimentally the case, the pump Rabi frequency f2 is con- 
stant along the vapor cell and these equations are simply 
first order coupled linear differential equations. When 
the dynamics is dominated by the cross-terms and no 
conjugate field is injected, the probe and conjugate fields 
grow asymptotically exponentially and the system be- 
haves like a phase-insensitive amplifier for a probe input 
field |T31 [TB] . At the quantum level, these cross-terms are 
responsible for the creation of quantum correlations be- 
tween the output probe and conjugate fields, leading to 
the production of a two-mode squeezed state [TJ I15j . The 
larger the gain of the amplifier, the greater the amount 
of squeezing. 

The full dynamics is more complicated than pure 4WM 
because of the presence of the direct terms Xpp an d Xcc- 
However the general form of Eqs. ^ and Q offers us 
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a straightforward interpretation of Xpp an d Xcc in terms 
of effective linear susceptibilities for the probe and the 
conjugate fields. Note that unlike the usual linear sus- 
ceptibilities, these effective susceptibilities depend, non- 
linear ly, on the pump field. Therefore they give rise to 
a pump-dependent complex index of refraction for the 
probe and the conjugate. The real part influences the 
phase-matching of the process, as discussed below. The 
imaginary part translates into absorption. 

In the original proposal by Lukin et al [5J, the pump 
beams are resonant with an atomic transition (Ai = 
A-2 = 0). This results in a remarkable situation where 
the pumps and the twin beams fulfill the two-photon Ra- 
man resonance and enter an electromagnetically-induced 
transparency (EIT) condition [17]. As a result, they see 
a perfectly transparent medium, as witnessed by a van- 
ishing imaginary part of Xpp an d Xcc- At the same time, 
the cross-susceptibilities are enhanced by the coherence 
of the hyperfine electronic ground states. In theory, it 
should lead to very efficient 4WM and virtually no ab- 
sorption, even for a weak pump. In practice, the EIT 
effect in hot vapors is limited because the Doppler effect, 
the presence of multiple excited hyperfine levels, and the 
finite transit time of the atoms in the laser beams all 
act to increase the decoherence rate between the ground 
states. This causes residual absorption of the probe field 
and as a result, only low levels of squeezing have been 
observed in this configuration |18j . 

In contrast, most recent squeezing experiments in hot 
atomic vapors operate at large detuning Ax/2tt, typically 
0.5 to 1 GHz, and at larger pump power [glfTOl E]. In 
these conditions, the susceptibilities take a different form 
from the resonant case. This off-resonant form is depicted 
in Fig. [2j for typical experimental parameters. 

The first important property of these susceptibilities 
is that the 4WM resonance is shifted from the bare two- 
photon resonance (6 = 0) by the light shift created by 
the more resonant pump, on the 5Si^(F = 2) — > 5Pi/2 
transition. For typical experimental parameters, the shift 
is of the order of —57. In the rest of the paper, we call 
"blue side" of the 4WM resonance the range of detunings 
6 for which the frequency of the probe is above the 4WM 
resonance. The other side of the 4WM resonance is the 
"red side" (see Fig.[2|. 

The second property is that the imaginary part of \pp 
(dashed line) is maximum at the 4WM resonance, due to 
Raman absorption. It is therefore not a good place to ob- 
serve 4WM because the medium is essentially opaque for 
the probe at this detuning. On the other hand, Q(Xpp) 
decays much faster than the magnitude of Xpc when mov- 
ing away from resonance, therefore there is a range of (5 
on each side of the 4WM resonance where Xpc still ex- 
hibits a substantial magnitude while Q(xpp) has almost 
completely vanished (see Fig. [3]). These regions of the 
two-photon detuning are better places to observe quan- 
tum effects. 

The third important property is the behavior of the 
real part of Xpp (solid line), which is effectively respon- 
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Figure 2. (color online) The direct and cross susceptibilities 
for the probe and conjugate fields as a function of the two- 
photon detuning S, varied by changing the probe frequency, 
expressed in units of the excited state decay rate 7 = 2tt x 6 
MHz. The solid black lines are the real parts, the dashed 
red lines are the imaginary parts. The units on the y axes 
are arbitrary but identical for all four susceptibilities. The 
resonant Rabi frequency of the pump is f2 = 6O7, the detuning 
of the pump is Ai = 1407, and the decoherence rate of the 
excited state is estimated (see below) to be y c = 0.2-y. The 
feature at S — I5O7, visible on Xvv i n the insert, corresponds 
to the one-photon resonance of the probe from the F — 2 
hyperfine ground state. The region of positive 5R(Xpp) around 
the 4WM resonance is indicated with a thicker line. 
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sible for the index of refraction for the probe. Around 
the 4WM resonance, $t(Xpp) is the sum of a disper- 
sive feature resulting from the 4WM coupling and the 
off-resonance negative susceptibility resulting from the 
one-photon transition between ground state and excited 
state. These competing terms lead to a cancellation of 
•ft(Xpp) around the bare two-photon resonance (5 = 0). 
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Figure 3. 4WM coupling jxpcj (solid black line) and probe 
absorption Ss(Xpp) (dashed red line) as a function of the two 
photon detuning S, for 7 C = O.27. The probe absorption is 
also shown for 7 C = 7 (dashed-dotted line) and for y c = O.O27 
(dotted line). In the latter case, the transparency window 
around S = is visible. 

It is legitimate to wonder if it would not be possible 
to take advantage of the coherence between the ground 
states to reduce the probe absorption through the EIT 
phenomenon, as envisioned in the original proposal. The- 
oretically, for a very low value of the ground state de- 
coherence rate 7 C , it is indeed possible to observe the 
transparency window at the bare two-photon resonance, 
as shown in Fig. [3j In practice, such low decoherence 
rates are not achieved due to imperfections associated 
with hot vapors. Experimental results reported below 
and in Refs [TH] and [50] are compatible with a higher value 
7 C ~ O.27, for which there is no marked transparency 
window. In spite of the lack of efficient EIT, the off- 
resonant value of 3(xpp) is small enough compared to 
the cross-coupling to ensure efficient 4WM, and exper- 
imental set-ups do not require special precautions with 
regards to decoherence sources such as stray magnetic 
fields. 

Note that these considerations do not apply to the con- 
jugate field, which is much further detuned from reso- 
nance than the probe. The direct susceptibility of the 
conjugate Xcc is substantially smaller than the other sus- 
ceptibilities and can be regarded as zero in practice. 

Now that adequate ranges of 8 have been identified 
for which probe losses are negligible compared to 4WM 
amplification, the question is whether there is a geometric 
configuration of the light fields for which 4WM is phase 



matched. To answer this question, we recall from Ref. 6 
the solutions to the propagation equations ^ and Q, 
for a seed £ s on the input probe and no input conjugate: 



£ p = £ s exp(SaL) 



cosh(££) + - sinh(£L) 



£* = £ s exp(<5a£)-^ smh(£i) 



(5) 
(6) 



as I £p ' 



where L is the length of the medium, a P j = ik p x P j/2, 
a c j = ik c x*j/2, Sa = (a pp ~ a cc + iAk z )/2, a = (a pp + 
a cc - iAk z )/2, £ = yj-a pc a cp + a 2 , and j = p,c. We 
define the probe and conjugate intensity gains g p and g c 
p \£ s \ 2 and \£ c \ 2 =g c \£ s \ 2 . 
From these expressions, we plot in Fig. [4] the probe 
and conjugate gains as a function of S and the geometric 
phase mismatch Ak z . One can see that for Ak z ~ 0, 
the maximum gain is obtained around <5 ~ 0, for both 
the probe and the conjugate. When Ak z increases, the 
gain on both the probe and the conjugate increases while 
the gain resonance moves towards the 4WM resonance 
(shown with the vertical dashed line). At larger Ak z , 
the gain resonance comes asymptotically within 7 of the 
4WM resonance, the probe intensity drops while the con- 
jugate intensity keeps increasing. Finally, at large Ak z , 
the conjugate intensity also drops. 
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Figure 4. Theoretical output probe and conjugate gains g p 
and g c as a function of the two-photon detuning and the 
geometrical phase mismatch. Here the decoherence rate is 
7c = O.57, the atom density is Af = 3 x 10 12 cm" 3 , the length 
of the medium is L = 12.5 mm and the pump Rabi frequency 
is Q = 6O7. The dashed lines indicate the position of the 
4WM resonance. 

The position of the gain resonance is the result of the 
4WM phase-matching (or absence thereof), and is in- 
fluenced by the effective index of refraction seen by the 
probe and the conjugate as follows. Since 5t(xcc) is much 
smaller than the other susceptibilities, the conjugate ef- 
fectively propagates in a medium of index f . The situa- 
tion is different for the probe, for which the index of re- 
fraction changes sign at the 4WM resonance and at 6 ~ 0, 
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as indicated in Fig. [2[ The change in sign of 5i(Xpp) at 
5 ~ means that the probe experiences an effective index 
of refraction, n p , smaller than 1 for 6 > and larger than 
1 between S < and the 4WM resonance. We assume 
that the index of refraction experienced by the pump, 
no j is close to 1 since the pump tends to optically pump 
the atoms towards the ground state of the off-resonant 
transition 5S 1 / 2 (^ = 3) — > 5Pi/ 2 - This assumption will 
be refined later. 

The geometric phase-matching condition Ak z = is 
the phase-matching condition in free space. It is fulfilled 
only when the beams are rigorously co-propagating, as 
shown in Fig. [5]ja). For the process to be efficient, the 
effective phase-matching condition must be fulfilled: 



2k n - 



ipkp k c — 0, 



(7) 



as shown in Fig. [5jb). This condition is identical to 
the geometric phase-matching condition (Afc z = 0) only 
when n p = v/T + $i(Xpp) — 1) which occurs around 6 = 
(Fig. [2]). When n p > 1, the effective phase-matching con- 
dition ([7]) imposes Ak z > 0, which corresponds to hav- 
ing a finite angle 9 between the pump and the probe and 
conjugate [21] . This occurs on the red side of the 4WM 
resonance, for 5 < (see Fig. [2 1. As 9 increases, the 
gain resonance is shifted towards higher values of 5ft(x pp ) 
and therefore towards the 4WM resonance. A negative 
geometric phase mismatch Afc z cannot be fulfilled. For 
this reason, no effective phase matching can happen on 
the blue side of the 4WM resonance, where n„ < 1. 



III. EXPERIMENTAL VERIFICATION 

In order to verify the theoretical predictions, a test 
was performed, as shown in Fig. [6] in which a 750 mW 
pump laser of beam waist 0.9 mm drives the Dl line 
at 795 nm in a 12-mm-long cell of 85 Rb vapor, heated 
and temperature stabilized at ~ 110°C. A seed beam at 
the probe frequency is produced by diverting a fraction 
of the pump through an AOM operating at cjhf/2 — 
2ir x 1.5 GHz in a double-pass arrangement. This seed 
beam, of power 10 — 20 /xW and waist 0.4 mm, then 
intersects with the pump inside the cell at a small angle 
9. The probe and conjugate beams are perpendicularly 
polarized with respect to the pump which is rejected at 
the output with a polarizing beam splitter. From the 
measured input seed power P s , output probe power P p , 
and output conjugate power P c , the probe and conjugate 
gains {g p = P p /P s and g c = P c /P s , respectively) are 
obtained. 
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Figure 5. (a): configuration where the geometric phase- 
matching condition is fulfilled (Ak z = 0); (b): configuration 
where the effective phase-matching condition for an effective 
index of refraction of the probe n p > 1 is fulfilled. In this case, 
there is a necessary geometric phase mismatch (Ak z > 0). 
The wave- vectors in vacuum ko, k p and k c , for the pump, 
probe and conjugate respectively, all have nearly the same 
magnitude. Energy conservation ensures that k p + k c = 2ko. 

Close to the 4WM resonance, the increase in ^s(xpp) 
accounts for the reduction in probe power with respect 
to the conjugate power, seen in Fig. |4j In this region, a 
high level of probe absorption coupled to a large 4WM 
gain still produces a strong conjugate output. At larger- 
angle 9, the effective phase-matching condition requires 
a value of 6 so close to the 4WM resonance that the high 
probe loss prevents the 4WM from happening at all. 



Figure 6. The Rb cell is pumped with a bright pump beam, 
in blue, and seeded at an angle 9 with a probe of power P a , 
in red. Emitted is an amplified probe of power P p , also red, 
and a conjugate of power P c , in yellow. The three powers, 
P s , P P and P c , are monitored by photodiodes while the probe 
frequency is scanned. PBS: polarizing beamsplitter. 

A dichroic-atomic-vapor laser lock [55] is in place to 
regulate the pump frequency thus maintaining a constant 
Ai. The angle of intersection 9 between the pump and 
probe is set by manually adjusting a pair of input mirrors, 
and 6 is adjusted by changing the AOM drive frequency. 
For a selection of values of 9, S/2n is swept over a range 
of typically 60 MHz, where noticeable gains g p and g c 
are observed. The angular range used, from 0° to 1°, 
satisfies the condition that there must be a full overlap 
of the beams over the cell length (this is achievable for 9 
up to 5°). Figures [7]ja, b) show the contour plots of g p 
and g c as a function of S and 9. 

The general features of figure [4] are reproduced, in- 
cluding the shift of the gain peaks towards the 4WM 
resonance when 9 is increased, as well as the crossover be- 
tween probe power and conjugate power. The crossover 
is the value of 9 for which the peak conjugate power is 
equal to the peak probe power. The main discrepancy be- 
tween the experimental data and the theoretical predic- 
tion is the measured drop in probe and conjugate power 
for 9 > 0.6°. This leads to peak gains for both the probe 
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and the conjugate which are well below the theoretical 
prediction. There are two main reasons for this behav- 
ior. Firstly, for a finite 9, the Doppler effect due to the 
thermal motion of the atoms does not cancel between 
the pump and the twin beams. For 6 = 1°, the resid- 
ual Doppler effect on S reaches 2j, which is roughly the 
width of the gain peak itself. This results in a broadening 
of the gain resonance at larger angles 9, when compared 
to the resonance given by the theoretical model, which is 
visible in Fig. [7J 

Secondly, and perhaps more importantly, when 9 is in- 
creased past 0.5°, the probe beam is subject to a strong 
effective cross-Kerr interaction with the pump around the 
gain resonance. This is due to the fact that as the gain 
peak moves closer to the 4WM resonance the effective 
index of the probe is resonantly enhanced, as shown by 
the behavior of 3?(Xpp) m Fig. [2j As a result, the trans- 
verse intensity variation of the pump realizes a strong 
lens for the probe, causing it to emerge from the cell 
with an angle of divergence comparable to or larger than 
9 itself [23]. 

The good agreement between the theoretical model 
and the measurements at angles where the probe fo- 
cusing is negligible gives us the opportunity to extract 
the values of those parameters which are not otherwise 
easily obtained, in particular the pump Rabi frequency 
fl, the atom density Af, and most importantly the de- 
coherence rate j c . Moreover, as explained in the ap- 
pendix, we can also introduce the effect of the index 
of refraction of the pump on the phase-matching con- 
dition by replacing the geometric phase mismatch with 
Ak z — 2nofc — k p cos 9 — k c cos 9, where no = 1 — £ is the 
pump index of refraction. By fitting the model on the 
data, as shown in Figs. [7^c, d), we find that fl = 6O7, 
N = 2.8 x 10 12 cm" 3 , 7c = O.27, and e = 6.5 x 1CT 6 . For 
this data set, the pump detuning was determined to be 
Ai = 1407 by calibrating the position of the 4WM gain 
against a Rb spectroscopy spectrum. As expected, the 
index of refraction for the pump on the blue side of the 
atomic resonance is smaller than one. 

One can check that the parameters extracted from the 
fit are broadly consistent with the estimated experimen- 
tal conditions. Our pump beam parameters lead to a 
peak intensity of 60 W • cm~ 2 , which results in a reso- 
nant Rabi frequency fl — 8O7 for a mean electric dipolc 
d = 1.47 x 10~ 29 C • m [23] . From the vapor pressure data 
summarized in Ref. [24j the number density of 85 Rb at 
100° C is J\f ~ 4 x 10 12 cm" 3 . Finally the index of refrac- 
tion for the pump, evaluated for ground state populations 
of 6% in the lower hyperfine state and 94% in the upper 
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hyperfine state, as given by Eqs. (A. 17 A. 20), and for the 
electric dipole value given above, is n a = 1 — 1.6 x 10~ 5 . 
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Figure 7. Probe gain g p and conjugate gain g c as a function 
of the two-photon detuning and the probe-pump angle. Top: 
experimental data. Bottom: fit to the experimental data for 
the angles between 0.1° and 0.5°, that is to say for the region 
below the dashed line. The hatched regions represent values 
outside the color scale. The theoretical peak gains for the 
probe and the conjugate are 350 and 1100 respectively. 



the atoms in the cell. The good agreement between the 
model and the experimental data suggests that consid- 
ering only average values of the single-photon detunings 
Ai and A2 captures most of the physics at play. How- 
ever, considering that the probe is typically tuned to the 
edge of the Doppler profile, it is legitimate to wonder 
what level of absorption this causes. It turns out that 
although the EIT has a limited impact, the pump field is 
highly saturating at resonance and causes a wide Autler- 
Townes splitting for those atoms resonant with the probe 
field. This renders even the resonant part of the atomic 
vapor highly transparent for the probe as long as the 
optical depth is not too large. In practice, noticeable 
levels of squeezing can be observed for Ai/27r as low as 
500 MHz [7], which is well inside the Doppler profile at 
our operational temperature. 

In order to verify these assumptions, we extended the 
above model to include the full Doppler distribution of 
detunings. The result, fitted to the gain curve of the 
probe as a function of S, is shown in Fig. [8] The small 
discrepancy in the width of the single-photon resonance 
is due to the fact that the model does not include the 
hyperfine structure of the excited state, whose main effect 
is to broaden the apparent Doppler profile. 



IV. DOPPLER BROADENING 



OPTIMIZING FOR QUANTUM NOISE 
REDUCTION 



The model developed above does not take into account 
the Doppler broadening caused by the thermal motion of 



The observation of large levels of intensity-difference 
quantum noise reduction requires a near-perfect phase- 
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Figure 8. Probe gain as a function of the two-photon detun- 
ing at low angle 6 and fit to the data including the Doppler 
broadening in the model. The discrepancy in the width of the 
Doppler-broadened single-photon absorption dip is caused by 
the presence of 2 excited hyperfine levels separated by 6O7. 
The model only considers two excited states which are degen- 
erate in energy (Fig. Si). 



insensitive amplifier with a gain of at least a few units. 
In the case of 4WM in a hot atomic vapor, this means 
that the absorption of the twin beams must be kept to 
a minimum while ensuring that efficient 4WM can take 
place. As pointed out previously, the two-photon de- 
tuning which fulfills those two conditions is S ~ 0, and 
not the 4WM resonance in itself. This is firstly because 
away from the 4WM resonance the probe susceptibility 
responsible for the absorption drops faster than the cross- 
susceptibility responsible for the 4WM gain, and secondly 
because the ground state coherence reduces the probe ab- 
sorption at that detuning [25] . 

From the previous discussion, at a detuning 5 ~ 0, 
effective phase matching of the 4WM process requires 
geometric phase-matching Ak z = 0. When corrected for 
the index of refraction for the pump, this condition cor- 
responds to the introduction of a small angle 9 -C 1° 
between the pump and the probe beams. This is indeed 
how the best levels of squeezing have been experimen- 
tally observed [7HIu]. The production of squeezed light 
by non-degenerate 4WM in a hot vapor therefore ben- 
efits from two favorable circumstances. First, the opti- 
mum angle 9 is small enough that the residual Doppler 
effect acting on the nearly co-propagating beams is much 
smaller than the width of the gain peak. Second, the 
same angle is large enough that the beams participating 
in the 4WM can be spatially separated at the output of 
the vapor cell. 

It is worth noting that in certain conditions, a small 
amount of loss on the probe beam can be beneficial. For 
instance, for a probe seed containing a large amount of 
classical noise, it is useful to have twin beams of equal 
powers in order to ensure proper rejection of the classical 



noise in the balanced detection 7J 26 . In particular, the 
existence of points in the (6, 9) parameter space where the 
probe and conjugate powers are perfectly balanced has al- 
lowed the detection of 8 dB of intensity-difference squeez- 
ing at frequencies as low as 2.5 kHz, despite the presence 
of substantial technical noise on the input probe [TJ. 



CONCLUSION 

We have shown that it is possible to phase-match 
4WM in a hot atomic vapor so that absorption is re- 
duced to a level where quantum effects could manifest 
themselves. Most reports of large squeezing generated 
by non-degenerate forward 4WM in hot vapors to date 
have indeed used a similar arrangement as the one pre- 
sented here, with nearly identical atomic and beam pa- 
rameters. Furthermore, the model we have developed 
appears very accurate for small angles between beams, 
even in its simple form neglecting both the Doppler ef- 
fect and the Zeeman sub-structure. It allows the robust 
extraction of parameters which would be difficult to de- 
termine by direct measurement, such as the decoherence 
rate of the ground states. 

By providing further insight into the mechanism of 
the process, the present results may lead to the real- 
ization of different configurations or different regimes 
amenable to the production of interesting quantum states 
of light, such as multi-spatial-mode phase-sensitive am- 
plifiers [571. 
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Appendix: Expressions for the susceptibilities 

1. Derivation of the susceptibilities for the probe 
and the conjugate 

In this appendix we derive the dynamics of the double- 
lambda configuration described in the paper. The re- 
sponse of an atomic system to an optical field is deter- 
mined by the polarization of the medium, which acts as 
the driving term in the wave equation. For a medium 
that consists of non-interacting particles, such as a di- 
lute atomic vapor, the polarization of the medium is of 
the form [13] 

P=AA(d), 

where Af is the number density of the atomic medium 
and d is the atomic dipole moment operator. The po- 
larization of the medium can be written in terms of the 



atomic eigenstates, such that it takes the form 



where the sum is over all the involved atomic transitions, 
^f,nm is the frequency of the field that couples the tran- 
sition between levels n and to, a nm is the density matrix 
element between levels n and to in a rotating frame at 
frequency u>f tTlm , and d„ m is the dipole matrix element 
between levels n and to. For an isotropic medium, the 



polarization of the atomic medium at a particular fre- 
quency is given by 

P(^/,nm) A/"d mn (7 nm . 

Thus, the response is completely determined by the 
atomic coherence of the corresponding transition. 

The equations of motion for the density matrix ele- 
ments in the rotating frame can be shown to be of the 
form [T3] 



a —i(A )rr + - \d ■ E(r f\a p- i ( t "/.<"»-'"/.t»«0* — n H ■ E(r t\ e - i 'y l »f-™-^f,™') t 

u nm — 1 \ L -*nm Jnm ) u nm ~ ^_ / ^ |^ u ni/ J -' V ) L ) u vra^ u nv^vra J -' V ; b j c 

v 

&nn — ~T ^ ^ [d n ,y ■ E(r, i)(7 l y n 6 ^ T nv^-vn ' E(l*, £)e ^ j ~t- ^ ^ F nm (7 mm ^ ^ ^mn&nn: 



E m >E„ 



E m <E n 



for n ^= to 



where A nm is the detuning of the field at frequency w/; nm 
from the transition between levels n and to , T mn is the 
population decay rate from level n to level to, 7„ m = 
(r„ + r m )/2 + 7^ m is the dipole dephasing rate, T n is 
the total decay rate out of level n, and 7^ m is the dipole 
dephasing rate due to any other source of decoherence. 

We now specialize to our 4WM process in the double- 
lambda configuration, shown in Figs, [I] and [9j with a 
single pump field, £o- For this case the total electric field 
is of the form 



E(r,t) = f e i(ko ' r "" Jot) e + £ c e l(kc 



where So, £ p , and £ c are the field amplitudes for the 
pump, the probe, and the conjugate, respectively, and 
Ei are unit vectors describing the polarization of the 
fields. We assume that the pump couples the transitions 
|1) — » 1 3) and |2) — > |4), the probe couples transition 
|2) — > |3), the conjugate couples transition |1) — > |4), 
and that the transitions |1) — > |2) and |3) — » |4) are not 
dipole allowed. With these assumptions and using the 
rotating- wave approximation, the equations of motion for 
the density matrix elements take the form 



-E p e i ^ r - U ^e p + c.c. 



& n = 




+ o; e - 
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where fli and O2 give the Rabi frequencies for the two 
transitions coupled by the single pump, tt p the Rabi fre- 
quency for the transition coupled by the probe, and Q c 
the Rabi frequency for the transition coupled by the con- 
jugate. In order to completely eliminate the explicit time 
dependence when doing the rotating-wave approxima- 
tion, the frequencies of the pump, probe, and conjugate 
fields need to satisfy the relation 

= (dp ) 

which is just an energy conservation condition for the 
4WM process. 



in steady-state condition for 041 and 032 to find that 




Figure 9. Double-lambda scheme with a single pump field 
Eq. States 1 3} and |4) are orthogonal linear combinations of 
magnetic states of the excited hyperfine levels. 

In order to solve the above equations and obtain ana- 
lytical expressions, we assume that the probe and conju- 
gate fields are weak fields, such that we only keep terms 
to first order in Q p and f2 c . In this case, the polarization 
of the medium at frequency u>i (where i indicates probe 
or conjugate frequency) can be divided into two different 
terms, one that is proportional to the field at frequency 
Wi and one that is proportional the field at frequency 
2wq — LUi, such that 



P(u p ) = e XppH)£ P e 4Vr 
P(uj c ) = eoXcc(w c )f c e 4kc ' r 



(A.11) 
(A.12) 



In doing this, we have introduced the susceptibility of 
the atomic medium Xiji which completely characterizes 
the response of the atomic system for a given field. The 
direct susceptibilities Xpp.cc act as the effective linear sus- 
ceptibilities for the probe and conjugate, respectively; the 
cross-susceptibilities Xpc,c P are responsible for the 4WM 
process. 

In addition to the approximations mentioned above, we 
assume the dipole moments for the two pump transitions 
to be equal, such that Oi = O2 = fi- Under these ap- 
proximations we can solve the density matrix equations, 
Eqs. ( A.lflA.lO ), to all orders in the pump field (Q) and 



Xpp 
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"022,44 



Xcp = 



^14C?23^32^ 2 \C 



(A.15) 



e hD\n\ 2 
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where we have defined 

= (e«+6i)(&4 

the population differences 

CT 11,33 = fll — 033 = 
011,44 = fH — 044 = 
022,33 = °22 — 033 = 
022,44 = 022 ^ 044 = 

and the complex decay rates 
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(A.18) 
(A.19) 
, (A.20) 
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= *(A 2 - 
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£32 
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= iAi - 


61 


= iS - 7c 



7 
2 
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In deriving these equations, we have assumed that the 
total decay rate out of the excited states are the same and 
with equal branching ratios to the two ground states; that 
is, r 4 = T 3 = 7 and V u = T 2 4 = T 13 = T u = 7/2. In 
addition, we assume that the additional dipole dephasing 
term 7?- is only significant for the ground state coherence, 
such that 7j C j^i2 = an d 7i2 = 7c- 
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2. Index of refraction for the pump 

In the above derivation of the susceptibilities for the 
probe and the conjugate, we have implicitly set £q, the 
amplitude for the pump electric field, to be a constant 
throughout the medium. This assumes that the pump is 
neither dephased or absorbed. In practice, the pump de- 
tuning is large enough to neglect absorption, but even a 
small index of refraction may have a substantial impact 
on the phase-matching condition of the 4WM process. 
The effect of refraction on the pump is to multiply So 
(and therefore O) by a running phase factor as the pump 
propagates. This can be taken into account by simply 



replacing ko in Ak by noko, where no is the index of 
refraction. The index of refraction created by the pop- 
ulation Mi of the ground state i can be estimated to be 
n o = VI + X — 1 + x/2 with 

N,d 2 A, 
X e ftA2+ 7 2/4' 

where d is the dipole matrix element of the transition for 
large detunings and A, is the detuning, which is taken to 
be much larger than the hyperfine splitting of the excited 
state. 
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